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Abstract 

Spectral density of current fluctuations in a short ballistic superconducting 
quantum point contact is calculated for arbitrary bias voltages V. Contrary 
to a common opinion that the supercurrent flow in Josephson junctions is 
coherent process with no fluctuations, we find extremely large current noise 
that is caused by the supercurrent coherence. An unusual feature of the 
noise, besides its magnitude, is its voltage dependence: the noise decreases 
with increasing V, despite the fact that the dc current grows steadily with V. 
At finite voltages the noise can be qualitatively understood as the shot noise 
of the large charge quanta of magnitude 2A/V equal to the charge transferred 
during one period of Josephson oscillations. 

PACS numbers: 74.50.+r, 74.80.Fp, 73.20 Dx 
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It is well established that the supercurrent in classical Josephson tunnel junctions does 
not fluctuate - see, e.g., 0. This result is consistent with the notion that the supercurrent is 
the coherent dissipationless flow of Cooper pairs, and therefore is quite appealing intuitively 
and has acquired the status of a widespread opinion. There seems to be no reason to suspect 
that the situation should be any different in mesoscopic contacts with a few transverse 
modes, since, for instance, the noise properties of contacts in the normal state are insensitive 
to the number of transverse modes. Both in classical P|| and quantum f3],||] normal point 
contacts the shot noise was predicted to be suppressed due to the Pauli-principle correlations 
between electrons, and such a suppression has been found in experiments ||[7|. This picture 
is modified only slightly in point contacts between normal metals and superconductors, or 
in fully superconducting point contacts at large voltages, where Andreev scattering leads 
to partial reflection of electrons in the contact region giving rise to a finite level of shot 
noise 0|§. The noise in this case is associated only with the excess current, and therefore 
represents only a small fraction of the full classical shot noise. 

The aim of this work is to suggest and prove by microscopic calculation that contrary 
to intuitive expectation, the supercurrent flow in quantum point contacts should generate 
noise which is extremely large on the scale of the classical shot noise. The noise arises from 
the interplay between quasiparticle scattering and the supercurrent coherence. The latter 
amplifies the randomness of the quasiparticle scattering by "attaching" to each quasiparticle 
a large charge which is coherently transferred through the point contact. 

We begin by considering a single-mode ballistic contact between two identical super- 
conductors. Such contacts can be formed, for instance, in superconductor/semiconductor 
heterosctructures and exhibit quantization of the supercurrent ||. The length d of the con- 
tact is assumed to be much smaller than the coherence length £, as well as the elastic and 
inelastic scattering lengths in the superconductors, so that all scattering inside the contact 
region can be neglected. We are interested in calculating the spectral density of current I{t) 
in the point contact: 
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where z is the coordinate in the direction of the current flow and z = corresponds to the 
position of the point contact. Taking t and t' in this definition to be on opposite branches 
of the Keldysh contour, we express the current correlation function 

Kj(h,h) = \(I(h)I(h) +/(fa)/(ti)> - (1(h)) (I(t 2 )) (1) 

as a product of two nonequilibrium Green's functions. In the quasiclassical approximation 
this expression reduces to the following ||: 

e 2 

Mh, h) = - ? E T *[9±(h, h)a z gi(t 2 , h)a z + g<(h, h)a z gl(h, h)a z ] , (2) 

where the Green's functions g are 2x2 matrices in the electron- hole space, Tr is taken over 
this space, and J2± is the sum over the two directions of propagation through the contact. 
Here and below, cr's denote Pauli matrices. 

From this point on, the calculation proceeds differently for vanishing and finite bias 
voltages. At vanishing voltage, the contact is in equilibrium and all Green's functions and 
current correlation function (||) depend only on the time difference r = h — h- in the 
frequency domain the equilibrium functions g are: 

gl(e) = (1 - f(e))p ± (e) , g<(e) = -f(e)p ± (e) , (3) 

where p(e) = gn(e) — gA( e ) is the matrix of the density of states, and /(e) is the equilibrium 
Fermi distribution. Fourier transforming the correlation function (|2|) and using eq. (|3]) we 
get the following expression for the spectral density of current fluctuations 

Si{u>) = ^J dre-^K(r) = £ j def(e)(l - f(e ± hu))Tr[p ± (e)a zP± (e ± hu)a z ] . 

(4) 

We shall see below that in order to get a meaningful result for Si(u) we need to keep 
a finite energy relaxation rate in the model. We find the density p(e) which includes a 
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finite relaxation rate generalizing the Kulik-Omelyanchuk theory |H|. We perform the 
calculation in real time. For a short constriction the retarded Green's function can be 
written as gf^ + gjR, where is the spatially-uniform part of the Green's function of the 
jth electrode: 



o (0) 

9jR 



I 

V 
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I = e + Z7x , A = A + 



and gjR is the non-uniform part which satisfies the equation ||IT| , |T2 

dg R , 



dz 



$m ,9r] ■ 



(5) 



Here 7^2 are inelastic scattering rates, S = (e 2 — A 2 ) 1 / 2 and is the Fermi velocity. The non- 
uniform part of the Green's functions should also satisfy two boundary conditions: g — > 
at z — > ±00 (inside the electrodes), and (71 + g± — g 2 + g^ at z — 0. 

From eq. @ we obtain that the solutions decaying at infinity have the following matrix 
forms: 

/ ■ + . \ / 1 ( ae ^)Ti\ 

(ae^ 1 )^ -1 

where a is the amplitude of Andreev reflection from the superconductor: 



9ir — A \ 
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\ae 



1 [ae 



V 



a(e) = (e-6)/A. 

Matching solutions in the two electrodes with the boundary condition at z = we get 
the total Green's function at this point: 
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9r,± 



S cos(<^/2) =f ie sin(9?/2) 



V 



I 



(6) 



e cos(y>/2) =f i<5 sin(yV2) , Ae^ 
-Ae - *** , -ecos(y3/2) ± i<5sin(^/2) 

where y?x = (<pi + <£>2)/2. The advanced function gA can be expressed in terms of gR, 
9a{.(Py) = — 9*r(— ¥?)■ Using this relation and eq. (Q), we see that although the Green's 
functions are not gauge-invariant (i.e., depend on </?s), the current spectral density (f|) is 
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gauge invariant. Therefore, to simplify the notations, we take below = without loss of 
generality. 

In the most interesting limit of small energy relaxation, 7^2 < A/S, there are two distinct 
regions of energy where the density of states p is non- vanishing: one is | e |> A, and the 
other is the vicinity of the two subgap states with energies ±eo = ±Acos((^/2). Expanding 
eq. (H) in small deviations from ±eo, we get that in the subgap region the density of states 



is: 



/x „^ r , m 7(eo)A sin(W2) 
p ± (e) = 2Re[ gR)± {e)}- 



( 1 ±A 



(e T e ) 2 + 7 2 M/4 



(7) 



V =Fl -I J 

Using this relation in eq. (JD, we get the subgap contribution to the current noise at low 
frequencies uj ~ 7: 

(1) = / J sin(yg/2) \ 2 7 (e ) 

' 1 J l v cosh(6 /2T) y ' 27r(^ 2 + 7 2 (eo))' 1 ' 

where J = eA/h, and 

/ \ 0( M C / \\ ^, 9(6^ -A 2 ) (6-6Q 3 COsh(6/2T) 

7(e) = 2( 7l (e) - - 72 (e)) = a / & ^ _ A . sin h((e - eO/2T) cosh(eV2T) " (9) 

Here a is a constant determined by the parameters of electron-phonon interaction. 

Equation (||) shows that there is a very large low-frequency noise associated with the 
supercurrent flow through a quantum point contact. Although the total noise intensity 
decreases with decreasing temperature, its zero-frequency density can actually increase at 
low temperatures because of the rapid decrease of 7. The noise has a very simple time- 
domain interpretation. Equation (|]) implies that the current which is considered to be a 
"dc" supercurrent does not flow continuously; rather it is a stochastic process with a typical 
realization shown in the inset in Fig. 1. Because of the quasiparticle exchange between the 
bulk electrodes and the two subgap states ±eo localized in the point contact, the system 
jumps between the state carrying positive current Iq sin(<^/2) and the state carrying negative 
current — Iosm(ip/2) with an average rate 7(eo). The jumps occur in such a way that the 
probabilities of finding a positive and negative current are, respectively, /(— eo) and /(eo). It 



is straightforward to check that the spectral density of such a process is indeed given by eq. 
(§). Thus, the noise (^j) is the two-level noise which is an inherent part of the supercurrent 
thermalization in the quantum point contacts. 

Equations and (|6|) also determine the noise due to the the states above the gap. As 
we will see below, this part of the noise has a much less singular behavior in 7 than the 
subgap noise Therefore in order to calculate it we can limit ourselves to 7 = 0. In this 
case we obtain from eqs. (|J) and @: 

S?\co) = 4^ / ^t 1 - tanh (^) tanh(^^)] M (e) M (e + e II e + ^ I + 

sgn(e)sgn( e + hu)A 2 cos 2 (^/2)) , „( c ) ee ^J^W^f^ ' (10) 

An interesting feature of the noise ( PD is that it is phase-dependent, despite the fact that 
the states above the gap do not contribute to the average supercurrent. For T ^> A and 
lu — 0, eq. ([U]) gives: 

l + cos(<£>/2) 



i)]. (id 



1 — cos(</?/2) 

The first term in this expression is the regular Nyquist noise as in the normal state, while 
the second term is the phase-dependent correction associated with the gap. Because of the 
usual BCS singularity in the density of states, the second term diverges weakly as <p — > 0. 
This divergence is removed by 7; at small but finite 7, the logarithm is limited by ln(A/7). 
At (p — 0, S} (cu) has a similar weak singularity as a function of lu, 

s i H = ^rl T - 7r ln — 



7T 



h l 2 w 



The last component of noise comes from the "interference" of the subgap states and 
states above the gap in eq. (|). In the limit 7 — > we get from this equation and eqs. (||) 
and (0): 



S? ] H = E e((e ± ^) 2 - A 2 )(( eo ± M 2 " A 2 ) 1/2 [l - tanh(|,) tanh(^^)] 

12) 



Comparison of eqs. ( [LTD and ( fP2|) with eq. (H) shows that the subgap noise dominates at low 
frequencies as long as the temperature is not too large on the scale of A, i.e., T -C A 2 

Total spectral density of current fluctuation Si(u) as a function of frequency calculated 
numerically from eqs. (^) and ([]) is shown in Fig. 1 for several values of the Josephson 
phase difference ip G [0, n]. The Si(u) for <p G [it, 2tt] can be found from the relation 
Si(u,ip) = Si(u, 2ir — ip) which follows from (f|) and @. Figure 1 shows that the subgap 
noise (||) indeed dominates the spectrum even at not too small 7. We also see that due to 
the interference term (|l^) the high-frequency threshold of Si(u) at low temperatures shifts 
down from 2A/H with increasing (p. In particular, at *p = tt, the noise at large frequencies 
starts at A/% in accordance with eq. (|T2|). 

Now we turn to finite voltages. In this case, all observable quantities like the current 
correlation function (|2]) are periodic in time t = (ti + t 2 )/2 with the period of the Josephson 
oscillations To = nh/eV. This implies that the Green's functions g should be periodic with 
period 2T , since their off-diagonal elements can change sign when shifted by T . Therefore, 
g can be expanded Fourier series: 



ikeVt/h 
k 



where t = t\ — t 2 - Averaging the correlator (g) over t and Fourier transforming it with 
respect to the time difference r as in the static case, we get the following expression for the 
spectral density of current fluctuations: 

2 

Si{u>) = -j^r £ / dtrt\gZ{e)a,g< k {e ± hu)a z ] . (13) 

k, ±uj 



In eq. (|T3"|) we have taken into account the fact that summation over the two directions 
of propagation in the correlator (fj) is equivalent, for a symmetric contact, to summation 
over positive and negative bias voltages ±V. Since the noise is an even function of V, this 
summations only contributes a factor of 2. 

From the known solution for Green's functions of a short superconducting point contact 



13| , [T^| and standard relations between the different Green's functions: 



g > = \{g K + g R -g A ), g K = \{g K - g R + g A ) , 

we obtain the Fourier components g^. 

gk = VkOz + Qkicr-y . (14) 
Here pk is non- vanishing for even k and: 

p>(e)=2(l-F(e-eV/2))A(e), p£(e) = -2F(e - eV/2)A(e) , k > , (15) 

k 

A{e) = [] a(e + leV - eV/2) , 
i=i 

while qk is non-vanishing for odd k and is given by the same eq. (|i~5|) with odd k. For negative 
fc, pfc = p*_ k , and gfc = Q-k- Function F has the meaning of non-equilibrium distribution of 
quasiparticles in the point contact: 

oo n 

F{e) = + E II I fl ( e - meV ) I 2 t/( e - ( n + !)eV0 - /(* - neV)] , (16) 

n=0 m=0 



Combining eqs. ( H ) and (|T^) with eq. ([03) we get a final expression for Sj(u) 



5 7 (cj) = — ^E / ^(e)(l-i 71 (e±M)[ 1 + 2ReEIl a ( e + /e ' l/ ) a *( e + /e ' 1/± M]- 

(17) 



Equation ([T7|) together with eqs. (Q) and (0) are the main technical results of our work. 
Combined, these equations give the spectral density of current fluctuations in a short ballistic 
constriction between two identical superconductors at arbitrary voltages. As we can expect 
from our calculations for V = 0, the most interesting limit at finite voltages is V <C A/e. 
In this case eq. (|17D can be simplified further. Expanding the amplitudes a(e) of Andreev 
reflection in small relaxation rates 7^2 and replacing the sums with the integrals in eqs. ([16]) 
and ( |P7D we obtain the spectral density of current fluctuations at low frequencies, u <C A/h: 

SM = £f a «(# - F(« ± M) jf /" ^|g= }x 



cos( — [arccos( ^ ) — arccos — ]) , (18) 

where the quasiparticle distribution function reduces to 

F(e) = f{e) - de' 9 4- exp{- f dv ^ } . (19) 
W JK ' J-A de' P1 Je> eWA 2 - v 2 

Equation ( |18"D is still a very general result which gives the spectral density of current fluc- 
tuations for arbitrary relation between the bias voltage, temperature, frequency, and energy 
relaxation rate. At large temperatures T ^> A the relaxation rate 7 (||) and distribution 
function F are constant for energies inside the gap, and the zero-frequency spectral density 
(0) is: 

Srm- eA ' ^(l + A 2 )+2 + 2e-^ ^ 

The noise (|20| ) has an unusual voltage dependence; it decreases monotonically with increasing 
voltage despite the fact that, as can be shown, the average current grows with voltage. 
Although the monotonic decrease of the noise intensity with the voltage is a characteristic 
feature of large temperatures T ^> hj, noise always decreases with voltage at V ^> hj/e 
(see inset in Fig. 2). Indeed, in this case, eq. Ql9|) shows that F(e) = F(A), and we get for 
the spectral density ([18]): 



c / x _ eA2 1 + cos(7rfiu;/eV) 

I{ ) 27r 2 ^cosh 2 (A/2T)V (1 - (Hoj/eV) 2 ) 2 ' 1 ' 

At small voltages V <C fr'j/e, and temperatures T < A, the noise is independent of the 
voltage V: 

5/(0) = • (22) 

n z h 7(0) 

Equations ([H]) and ( p2|) describe two sides of the noise peak with the maximum at V ~ hj/e. 
The exact shape of this peak depends on the energy dependence of the relaxation rate 7 
and is shown in the inset in Fig. 2 for several temperatures in the approximation of energy- 
independent 7. The curves were calculated numerically from eq. fllTl). The main part of 
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Fig. 2 shows the zero-frequency spectral density of current fluctuations at arbitrary voltages. 
It illustrates a transition from the noise peak at small voltages to the large voltage regime 
where 5/(0) saturates at 4e 2 A/157r 2 fr. 

The fact that the noise at small voltages is very large on the scale of the regular shot noise, 
together with an unusual voltage dependence, reflects an unusual physical mechanism of the 
noise. For V 3> hj/e this mechanism can be described qualitatively as follows. The giant 
noise arises since each quasiparticle getting into the constriction region with energy equal to 
one of the gap edges generates an avalanche of Andreev reflections before it can escape out 
of the constriction by climbing up or down in energy to the opposite edge of the energy gap. 
The number of generated Andreev reflections is 2A/eV, so that each quasiparticle causes 
a coherent transfer through the constriction of a charge quantum of magnitude 2A/V '. For 
small voltages V< A/e this is much larger than the charge of individual Cooper pairs. In 
this way, the randomness of the quasiparticle scattering (quasiparticles get inside the energy 
gap with probability /(—A) from one electrode and with probability /(A) from the opposite 
electrode) is amplified. Therefore, the noise described by eq. fl2"T| ) can be interpreted as the 
shot noise of these large charge quanta, and is in fact the consequence of the coherence of 
the supercurrent flow through the point contact. 

It is interesting to note that this picture in the energy domain has a "dual" formulation in 
the time domain in terms of the non-equilibrium occupation of the two subgap states which 
are responsible for the dc supercurrent. In particular, the avalanche of 2A/eV Andreev 
reflections triggered by a quasiparticle corresponds in the time domain to one period of 
Josephson oscillation, during which the supercurrent Jo sm (v 9 /2) carries the charge 2A/V 
through the point contact. More generally, all small-voltage results for the spectral density 
of current fluctuations obtained above (eqs. (p0|) - (|22]) for the ac regime and eq. (H) for 



the dc regime) can also be obtained from a purely classical rate equation [14]] for the non- 
equilibrium occupation probabilities of the two subgap states. An advantage of such simple 
classical approach is the possibility to generalize it straightforwardly to time- dependent 
voltages and arbitrary bias conditions of the point contact. Of course an advantage of the 
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fully microscopic approach used in this work is that it gives spectral density of current 
fluctuations for frequencies and bias voltages that are arbitrary on the scale of A. 

The simple time-domain interpretation of the supercurrent noise at small bias voltages 
allows us to propose a simple generalization of our results for ballistic junctions to junctions 
with arbitrary transmission coefficient D. Indeed, making a very natural assumption that 
we have the same exchange of quasiparticles between the two subgap states localized in the 
point contact and the bulk superconductors, we get immediately for V = (cf. eq. (||)): 



Here N is the number of propagating transverse modes in the junction, = 
(eA/2h)Dksm((p)/[l — D k sin 2 ((/?/2)] 1 / 2 coincides at T = with the supercurrent carried 
by the kth mode, and e k = A[l - D k sm 2 (<p/2)] 1 / 2 00 ■ 

Since fluctuations of the current 
in all modes are independent, the current noise Si(u), expressed in terms of the total super- 
current, decreases roughly as 1/N with increasing number N of transverse modes. Therefore 
the giant small- voltage noise discussed in our work disappears in the classical limit iV — > oo. 
However, it can be very important in quantum point contacts with few transverse modes, 
and even in regular tunnel Josephson junctions of ultrasmall area. 

In conclusion, we have shown that the supercurrent flow in quantum point contacts leads 
to a giant current noise both in the regime of dc and ac Josephson effects. The noise arises 
from the interplay of quasiparticle scattering and coherence of the supercurrent flow and has 
an unusual temperature and voltage dependence. 
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(23) 



11 



REFERENCES 

[1] D. Rogovin and D.J. Scalapino, Ann. of Phys. 86, 1 (1974). 

[2] I.O. Kulik and A.N. Omel'yanchuk, Sov. J. Low Temp. Phys. 10, 158 (1984). 

[3] V.A. Khlus, Sov. Phys. - JETP 66, 1243 (1987). 

[4] G.B. Lesovik, JETP Lett. 49 592 (1989). 

[5] M. Biittiker, Phys. Rev. Lett. 65, 2901 (1990). 

[6] Y.P. Li, D.C Tsui, J.J. Heremans, J. A. Simmons, and G.W. Weimann, Appl. Phys. 
Lett. 57, 774 (1990). 

[7] M. Reznikov, M. Heiblum, H. Shtrikman, and D. Mahalu, Phys. Rev. Lett. 75, 3340 
(1995). 



J. P. Hessling, V.S. Shumeiko, Yu.M. Galperin, and G. Wendin, |cond-mat 950901!!. 



[9] H. Takayanagi, T. Akazaki, and J. Nitta, Phys. Rev. Lett. 75, 3533 (1995). 
[10] I.O. Kulik and A.N. Omel'yanchuk, Sov. J. Low Temp. Phys. 3, 459 (1977) 
[11] A.I. Larkin and Yu.N. Ovchinnikov, Sov. Phys. - JETP 41, 960 (1976). 
[12] A. Zaitzev, Sov. Phys. - JETP 51, 111 (1980). 
[13] U. Gunsenheimer and A.D. Zaikin, Phys. Rev. B 50, 6317 (1994). 
[14] D. Averin and A. Bardas, Phys. Rev. B, to be published; |cond-mat 9509153 



[15] A. Furusaki and M. Tsukada, Physica B 165&166, 967 (1990). 
[16] C.W.J. Beenakker, Phys. Rev. Lett. 67, 3836 (1991). 



12 



FIGURES 



Figure 1. Spectral density of current fluctuations in a short single-mode constriction 
between two superconductors in the regime of the dc Josephson effect. From bottom to 
top, the curves correspond to <p = 0; tt/2; it. The peak at low frequencies is the subgap 
contribution (|8|). The inset shows a typical realization of the "dc" supercurrent as a function 
of time at non-vanishing temperatures. The sign of the current switches randomly with the 
characteristic rate 7 giving rise to the low-frequency peak in spectral density. 

Figure 2. Zero-frequency density of current fluctuations in the ac regime as a function of 
the bias voltage at zero temperature. The inset shows a blow-up of the small-voltage peak 
at (from bottom to top) T = 0; 0.1 A; 0.2 A. 
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